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£f~^ • We study in detail the threshold energy dependence of the seasonal variation 

effect in the energy integrated solar neutrino signal of the Super-Kamiokande 
Q\ ■ detector in the case of the v e <-» v^ T vacuum oscillation (VO) solution of the 

solar neutrino problem. We show, in particular, that for the values of Am 2 
and sin 2 26 from the VO solution region, the predicted time and threshold 
e~ energy (T ej xh) dependence of the event rate factorize to a high degree of 
accuracy. As a consequence, the VO generated seasonal variation asymme- 
try is given by the product of an time-independent function of T e ,Th and the 
^ ■ standard geometrical asymmetry. For any given Am 2 and sin 2 29 from the 

■ VO solution region there exists at least one value of T e xh from the interval 

(5 — 11) MeV, for which the seasonal variation effect in the solar neutrino 
sample of events, formed by recoil electrons with kinetic energy T e > T e xh> 
is either maximal or very close to the maximal; it can vary dramatically with 
T e ,xh- One can effectively search for the VO induced seasonal effect by form- 
ing a set of samples of solar neutrino events, corresponding to a sufficiently 
large number of different values of T ei xh from the indicated interval and by 
measuring the seasonal variation in each of these samples. Predictions for 
the magnitude of the seasonal effect in such samples are given for a large set 
of representative values of Am 2 and sin 2 29 from the VO solution region. 
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1. Introduction 



The hypothesis that the solar neutrinos undergo vacuum oscillations into muon and/or 
tau neutrinos, v e «-> f M)T , when they travel from the Sun to the Earth continues to be a 
viable and very appealing explanation |8||| of the observed deficit of solar neutrinos [] [p!0|-p!5[ . 
The mean event rate data from the solar neutrino experiments Homestake, Kamiokande, 
SAGE, GALLEX and Super-Kamiokande, can be described in terms of two-neutrino v e <-> 
VfiM oscillations if the values of the two parameters, Am 2 and sin 2 29, characterizing the 
oscillations, lie in the region PJ1: 

5 x Kn n eV 2 £ Am 2 £ 5 x 1(T 10 eV 2 , (1) 

0.6 £ sin 2 29 < 1.0. (2) 

Although this result is obtained utilizing the standard solar model predictions of ref. |16| 
for the fluxes of the pp, pep, 7 Be, 8 B and CNO neutrinos, it is rather stable with respect to 
variation of the fluxes within their estimated uncertainty ranges J7|||. If, for instance, one 
treats the total 8 B neutrino flux as a free parameter in the mean event rate data analysis, the 
solution regions located within the limits determined by eqs. (1) and (2) change somewhat 
their position and magnitude only. In this case two new "low 8 B neutrino flux" solution 
regions appear (at 95% C.L.) at sin 2 29 ^ 0.6 for Am 2 ~ (3 - 4) x 10~ n eV 2 and Am 2 ~ 
(4 - 8) x 10~ 12 eV 2 (for a recent analysis see §|). For values of Am 2 ~ (4 - 8) x 10~ 12 eV 2 
the effects of the vacuum v e <-> v^m oscillations on the flux of 8 B neutrinos with energies 
E u ^ 5 MeV, which is presently studied in the Super-Kamiokande experiment and will be 
studied in the near future with the SNO detector, are hardly observable. This low 8 B flux 
vacuum oscillation (VO) solution can be tested in the future solar neutrino experiments 
BOREXINO, HELLAZ, etc., and we will not discuss it further. 

When the current Super-Kamiokande results on the recoil-e _ spectrum || are added in 
the analysis of the solar neutrino data, the regions of values of Am 2 and sin 2 29 for which 
the vacuum oscillations provide a good quality fit of the data, are considerably reduced and 
the best fit is achieved for Am 2 ^ (4.3 - 4.4) x 10" 10 eV 2 and sin 2 29 ^ 0.9. Actually, the 
present Super-Kamiokande data on the spectrum of the recoil electrons from the solar ( 8 B) 
neutrino induced reaction v + — > v + e~ favors the u e «-> iv T ) vacuum oscillation solution 
of the solar neutrino problem over the MSW solutions H. 

A strong evidence (if not a proof) that the solar neutrinos take part in vacuum oscillations 
on the way to the Earth would be the observation of a seasonal variation effect which differs 
from the standard geometrical one (see the second and the third articles quoted in [|IJ as well 
as Q-[7|,[r7|-£np ) . For values of Am 2 from the VO solution region (1) and neutrino energies 
E„ ^ 5 MeV, the oscillation length in vacuum, L v = 47rE v /Am 2 , is of the order of, or exceeds 
but not by a very large factor, the seasonal change of the distance between the Sun and the 
Earth, AR = 2eR , where e = 0.0167 is the ellipticity of the Earth orbit around the Sun and 
Ro = 1.496 x 10 8 km is the mean Sun-Earth distance. This leads to a noticeable dependence 
of the neutrino oscillation probability on the Sun-Earth distance, which in turn can create a 
specific seasonal difference in the solar neutrino signals in detectors like Super-Kamiokande, 
SNO and ICARUS. This difference should be particularly large in the signals due to the 



lr The possibility of solar v e <-> v s oscillations, v a being a sterile neutrino, is disfavored by the existing solar 
neutrino data |||i],£§. i 



0.862 MeV 7 Be neutrinos [|], to be studied with the BOREXINO detector. No analogous 
effect is predicted in the case of the MSW solution of the solar neutrino problem [fffi] . 

The seasonal effect due to the vacuum oscillations of solar neutrinos can be amplified or 
reduced by the standard geometrical one caused by the 6.68% decrease of the solar neutrino 
flux at the Earth surface when the Earth moves from perihelion (taking place in January) 
to aphelion (reached in July) on its orbit around the Sun ||. In general, the VO induced 
seasonal effect in the solar neutrino signals in the Super-Kamiokande, SNO and ICARUS 
detectors is predicted to be relatively small for values of Am 2 from the interval (1). For 
the envisaged and already attained rather low threshold e~-kinetic energy of detection of the 
recoil electrons in the Super-Kamiokande experiment T e Th — 5 MeV ||, the relative seasonal 
change of the event rate caused by the vacuum oscillations alone cannot exceed (even under 
the most favorable experimental conditions) approximately 10% || (see further). This is 
below the sensitivity currently reached in the search for this effect in the Super-Kamiokande 
experiment. For the v e <-> v^( T ) oscillations of interest the effect in the Super-Kamiokande 
detector is reduced, in particular, by the weak neutral current interaction contribution of 
the v^t) to the solar neutrino induced event rate. Correspondingly, the VO induced seasonal 
variation of the solar neutrino signal should be, in principle, somewhat larger (by a factor ~ 
1.2) in the charged current event rates in the SNO and ICARUS experiments. Nevertheless, 
the increase is not dramatic and it is worthwhile considering possible strategies which might 
lead to the enhancement of the seasonal effect in the samples of solar neutrino events collected 
by these or other real time detectors studying the 8 B neutrino flux. 

It was noticed in ref. that the magnitude of the seasonal effect generated by the vacuum 
solar v e oscillations in the Super-Kamiokande solar neutrino signal, and more generally in 
the solar neutrino signals in detectors based on the reaction u + e~ — > z/ + e~, is very sensitive 
to the change of the threshold e~-kinetic energy from T e> xh — 5 MeV to T 0j Th — 7 MeV. The 
VO induced seasonal effect in the signals in the SNO and ICARUS detectors should exhibit 
similar strong dependence on the chosen minimal value of the energy of the detected final 
state e~. In the present article we explore this observation. We study in detail the threshold 
energy dependence of the seasonal variation effect in the sample of solar neutrino events in 
the Super-Kamiokande detector in the case of the vacuum oscillation solution of the solar 
neutrino problem. We show, in particular, that for the values of Am 2 and sin 2 29 from the 
VO solution region (1) - (2), the predicted time and threshold e~ kinetic energy (T 6j Th) 
dependence of the energy integrated solar neutrino induced event rate factorize. The VO 
generated seasonal variation asymmetry is given by a product of a function which depends 
on T e Th but does not depend on time, and of the standard geometrical asymmetry. We show 
also that for any given Am 2 and sin 2 29 from the VO solution region there exists at least one 
value of the threshold energy T 6] Th from the interval (5 — 11) MeV, for which the seasonal 
variation effect in the solar neutrino sample of events, formed by recoil electrons with energy 
T c > T Cj Th; is either maximal or very close to maximal. This suggests a possible strategy to 
search effectively for the VO induced seasonal effect: it consists in forming a set of samples 
of solar neutrino events, corresponding to a sufficiently large number of different values of 
T 6i Th from the above interval, say, to T Gj Th = 5; 6; 7;...; 11 MeV, and by measuring the 
seasonal variation in each of these samples. Predictions for the magnitude of the seasonal 
effect in each of the indicated samples for a large set of representative values of Am 2 and 
sin 2 29 from the VO solution region are also given. 
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2. The Seasonal Variation Effect Observables 



We shall consider the simplest case of solar v e <-> io T ) oscillations generated by two- 
neutrino mixing. The probability that a solar electron neutrino with energy E v will not 
change into Z/V r ) on its way to the Earth when v e <-> v^u-) oscillations take place, has the 
well-known form: 

P(i/ e ^x/ e ;R(t),E I ,)=P(R(t))=P(t) = 1- J sin 2 20 [1 - cos2t3S], (3) 

2 L v 

where L v = 47rE„/Am 2 is the oscillation length in vacuum, 

R(t) = Ro [1 - e cos 2tt^] = R R(t), (4) 

is the Sun-Earth distance at time t of the year, T = 365.24 days, Ro = 1.496 x 10 s km, 
R(t) and e = 0.0167 being the mean Sun-Earth distance, the Sun-Earth distance at time t 
in A.U. and the ellipticity of the Earth orbit around the Sun. 

It is not difficult to check that in the neutrino energy interval of interest, 
E v = (5.0 — 14.4) MeV, the probability P(t) has, as a function of E Vl i) one minimum 
for values of Am 2 = (0.5 — 0.8) x 10~ 10 eV 2 , ii) one maximum and one minimum for 
Am 2 S (0.8 - 2.0) x lO" 10 eV 2 , etc.; for Am 2 (4.3 - 4.4) x 10~ 10 eV 2 it has 4 max- 
ima and 3 minima. If, for instance, Am 2 = 0.75 x 10 -10 eV 2 , P(t) decreases mono- 
tonically for E u rising in the interval E u = (5.0 — 9.05) MeV, and increases with E u for 
E u = (9.05 — 14.4) MeV; for Am 2 = 2.0 x 10~ 10 eV 2 it decreases as E u increases in the two 
intervals E u = (5.0 - 8.04) MeV and E u = (12.06 - 14.4) MeV, and increases with E u in the 
interval E v = (8.04 — 12.06) MeV. As it is easy to see from eq. (3), the derivatives of P(t) 



with respect to E v and R(t) have opposite signs 1 18(1 . Correspondingly, P(t) will decrease 
from January to June in the neutrino energy intervals, in which P(t) increases with E u , and 
it will increase from January to June in the energy intervals where P(t) decreases with E u , 
the seasonal effect having opposite signs in the two cases. There is no seasonal change of P(t) 
in the points of the extrema. Since the magnitude of the seasonal variation effect due to the 
vacuum oscillations is determined by the seasonal change of the oscillation probability P(t), 
it should be clear from the above discussion that for the VO solution values of Am 2 from 
(1), the seasonal effect in the solar neutrino induced event rate integrated over the entire 
energy interval E u = (5.0 — 14.4) MeV can be reduced considerably as a result of the mutual 
compensation between the opposite sign seasonal effects in the E„ subintervals where P(t) 
increases and decreases with E u . As we shall see, for certain values of Am 2 the compensation 
in the signal of the Super-Kamiokande detector is practically complete. We shall demon- 
strate in what follows that the indicated partial or complete compensation can be avoided 
in at least one of the samples of events corresponding to different minimal values of the 
recoil-e _ kinetic energy, i.e., threshold energies T 6i Th, from the interval T e = (5 — 11) MeV. 

The solar neutrino induced event rate in the Super-Kamiokande detector at time t of the 
year can be written for fixed threshold energy T 6) Th, Am 2 and sin 2 29 in the form: 

R(t,T e , T h;Am 2 ,#) = R(t,T CiTh ) = < J^— [r„ + (1 - r„)P(t)] > , (5a) 

R 2 (t) 
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where by the "average" of a quantity X, < X >, we shall understand in what follows 



< X > = ^ J dT c J dE„ X n(E„) (da(z/ c e~)/dT c ). (56) 

° T c , T h T e (l+^) 

Here F(B)/Rq = $(B) is the mean annual 8 B neutrino flux at the Earth surface, n(E„) 

14.4 MeV 

is the normalized spectrum of 8 B neutrinos, / n(E v )dE I/ = 1, dcx(z/ie~)/dT c is 

o 

the differential cross-section of the process v\ + e~ — > Vi + e~, I = e,/i (r), and 
r u = (dcr(z/ At e~)/dT e )/(da"(i/ e e~)/dT c ). In the neutrino energy interval of interest one has 
t u = (0.155 — 0.160). Expression (5b) is valid for ideal e~ detection efficiency and energy 
resolution of the Super-Kamiokande detector 0. It can be trivially modified to include the 
latter. The time-dependent quantities in the expression for R(t, T e xh) are the Sun-Earth 
distance, i.e., the geometrical factor (expressed in A.U.), R~ 2 (t), and the probability P(t). 

In what follows we shall consider the event rate R(t;T e! Th), averaged over a time in- 
terval with a central point t = t c and width At, R(t c , At; T 6i Th)- The point t c can be the 
time at which the Earth reaches the perihelion, t c = t p = 0, or aphelion, t c = t a = T/2, or 
any other chosen time of the year. It is convenient to choose t p = < t c < t a = T/2 since 
R(t c ) = R(T — t c ). The seasonal effect due to VO should exhibit a similar symmetry. We 
shall present results for t c = t p ; t p + T/12; t p + 2T/12; ...; t a and a width of one month, 
At = T/12. In their analysis of the solar neutrino data the Super-Kamiokande collaboration 
is utilizing time bins having a width of 1.5 months, At = T/8 ||. As we shall see, the choice 
of At ^ T/4 tends to suppress the seasonal variation effect. 

Obviously, the time-averaged event rate R(t c , At; T e Th ) is determined by the time- 
averaged probability P(t). It proves useful to consider two different types of time-averaging. 
The first includes the effect of the geometrical factor, 

tc + iAt 
1 f Pftl 

P GF(tc ,A t ) = - / dtgii (6, 

t c -|At V ' 



while in the second the geometrical factor effect is essentially eliminated: 

5 GF(t c ,Al 
K(t c ,At) 



P NGF (t c ,At) = ^#l^, (7) 



where 



tc+^At dt 



, s 1 f dt / t c \ sm27rAt/ 2T , 

/c(t c ,At) = — / =l + 2ecos 2tt^ ./ ; + OU 1 

V ; At J R2( t ) V TJ 2vrAt/(2T) V 

tc-^At 



Obviously, expressions (5a) and (5b) and the results based on them we shall obtain, are valid for any other 
"ideal" detector utilizing the v + e~ — > v + e~ reaction for detection of the 8 B neutrinos. In this sense our 
results can serve as a guidance for the expected T Cj Th — dependence of the magnitude of the seasonal effect 
in any detector of this type. 
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and we have used eq. (4) and the fact that terms ~ e 2 = 2.8 x 10~ 4 are beyond the 
sensitivity of the operating and planned solar neutrino experiments and can be neglected. 
It follows from eqs. (4) and (6) - (8) that the one year averaged probabilities P GF (t c , T) and 
Pngf^cjT) practically coincide: 

PGF(t C ) T) = P NGF (t c ,T)+0(e 2 ). (9) 

The same result holds for the average probabilities independently of At for 
t c = T/4; 3T/4, i.e., at the spring and autumn equinoxes: PcF(t c = T/4 (3T/4),At) = 
P NGF (t c = T/4 (3T/4), At) + 0(e 2 ). 

The expression for the average event rate in the time interval [t c — At/2, t c + At/2], 
in which the effect of the geometrical factor is essentially excluded, Rngf^o At; T e( Th), 
can formally be obtained from eqs. (5a) - (5b) by replacing the probability P(t) with the 
probability Pngf^c, At) and by setting the factor R~ 2 (t) to 1: 

R NGF (t c ,At;T CiT h) = < r„ + (1 - r,)P NGF (t c , At) > . (10) 

It is easy to see from eqs. (5) - (8) and (10) and the above remark that for the average event 
rate containing the contribution of the geometrical factor, R GF (t c , At; T e Th ), one has: 

R GF (t c , At; T e , Th ) = /c(t c , At) R N GF(t c , At; T e>Th ). (11) 

The one year averaged event rates R GF (t c , T; T e Th ) and Rngf^c, T; T c Th ) coincide up to 
corrections ~ e 2 . For t c corresponding to the perihelion, t c = 0, we shall denote them both 
by R(T ejT h): R(T e ,Th) = Rngf(0, T; T ejT h) — Rgf(0, T; T e ,Th). 
We shall analyze in what follows two observables: the ratio 

NiCtc, At; T e , Th ) = R i(^;T Th ) . = NGF) (12) 

and the related seasonal variation asymmetry, 

Ap as (t c , At; T c , Th ) = Ni(t c , At; T e , Th ) - N ; (t c + T/2, At; T e , Th ) , i = NGF, GF. (13) 

Since Ni(t c + T/2, At; T Ei Th) = Ni(T/2 — t c , At; T C) Th), it is sufficient to consider 
0< t c < T/4 in the case of the asymmetry. Note that both observables (12) and (13) do 
not depend on the total flux of 8 B neutrinos. The observable Nngf^c, At; T eiF h) coincides 
with the one introduced in ref. 0. The seasonal variation asymmetries which have been 
analyzed in refs. [|TT|-|2"1"|] are analogous to, but differ somewhat from, A^ GF / GF s (t c , At; T Gj Th)- 
In the absence of neutrino oscillations P(t) = 1 and we have: NN GF (t c , At; T eiF h) = 
1, N GF (tc,At;Te,Th)=K(t c ,At), A^ G s F (t c ,At;T CiTh ) = 0, and A s G e f; s (t c , At; T CiTh ) = 
Ageom(tc, At), where Ag^^(t c , At) is the seasonal asymmetry of purely geometrical origin, 

Areo & m (tc,At) = K(t c ,At)-K(t c + T/2,At) = 4 e cos(2vr^) ^^^^ + 0(e 2 ). (14) 



Note, that the perihelion was reached on 2, 5, and 4 of January in 1997, 1998 and 1999, respectively. 
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For < t c < T/4 (and At < T) the geometrical asymmetry is positive: A|g^(t c , At) > 0. 
As can be shown, the following simple relations hold true up to corrections ~ 10~ 3 if the 

'Mr) 



solar v e take part in two-neutrino oscillations v e <-> VuM with Am 2 ^ 5 x 10 10 eV 2- 



N NGF (t c , At; T Cj Th) = 1 + \ W(T CiTh ) A^ m (t c , At) + 0(1(T 3 ), (15) 



N GF (t c , At; T e , Th ) = N NGF (t c , At; T c , Th ) + -A s g ^ m (t c , At) + 0(1(T 3 ) 

= 1H-| (1H- W(T e , Th )) A- a o S m (t c , At) + 0(1(T 3 ), (16) 

A^ s F (tc, At; T e , Th ) = W(T e , Th ) A^ s m (t c , At) + 0(1(T 3 ), (17) 

Ag|f (t c , At; T c , Th ) = A^ F (t c , At; T c , Th ) + A^ m (t c , At) + 0(1(T 3 ) 

= (1 + W(T e , Th )) A-- m (t C) At) + O(10- 3 ). (18) 

The factor W(T e Th ) in eqs. (15) - (18) does not depend on the time variables t c and At and 
is given by: 

1 . < (l-r^TrStt sin27r^ > 
W(T«, Tt ) = - ^ 20 , (19) 

where P(Ro) is the one year average probability, i.e., P(R(t)), eq. (3), in which the Sun- 
Earth distance R(t) is replaced with the mean Sun-Earth distance Ro- In the case of absence 
of oscillations we have sin 2d = and/or Am 2 = 0, and therefore W(T e ,Th) = 0. 

The relations (15) - (18) can be derived using the following observations. For Am 2 ^5x 
lCT 10 eV 2 and E u > 5 MeV we have 27reRo/L v ^ 0.63. Correspondingly, one can expand 
the probability P(R(t)) in power series of x(t) = 27reR /L v cos27rt/T ||. In this way one 
obtains: 

P(R(t)) = P(R )+AP(t), (20) 

where 



AP(t) = - sin 2 29 
v ; 2 



(x(t) - lx 3 (t) + ...) sin27r^ - i(x 2 (t) - lx 4 (t) + ...) cos 2^ 



(21) 

For Am 2 < lO" 10 eV 2 and E u > 5 MeV one finds x(t) < 0.13. Thus, the term linear in 
x(t) gives the dominant contribution in AP(t): for the contribution of the quadratic term, 
for instance, we get 0.25x 2 (t) < 4.3 x 10~ 3 . Keeping only the term linear in x(t) in the 
expression for AP(t) one arrives at the relations (15) - (19). The linear approximation, as 
numerical studies we have performed showed, turns out to be equally accurate for 10~ 10 eV 2 < 
Am 2 < 5 x 10" 10 eV 2 . The reason for this somewhat unexpected result lies in the fact 
that for Am 2 > 10 -10 eV 2 and E u ~ (5 — 14) MeV, the argument of the t— independent 
sine and cosine functions in AP(t) is relatively large: 27rR /L v > (2.7 — 7.5). Thus, for 
Am 2 > 10~ 10 eV 2 , the functions (27reR /L v ) n sin2vrR /L v and (27reR /L v ) n cos 2ttR /L v 
with n > 2 are fastly oscillating functions of E u and the integration over the neutrino energy 
renders them negligible: 

< -^-(27reR /L v ) n sin(27rR /L v ) > ~ O(10^ 3 ), n > 3, (22a) 
2 n! 
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< -^-(27reRo/L v ) n cos(27rR /L v ) > ~ O(10~ 3 ), n > 2. (226) 
2(n!) 

Relations (17) and (18) imply that to a high degree of accuracy the VO generated sea- 
sonal variation asymmetries A^qp GF (t c , At; T e) Th) are proportional to the geometrical one, 
Aj|gf(gf)(^c, At; T e ,Th) ~ Aj^ s m (t c , At). Actually, the time (t c and At) and the energy de- 
pendencies of the observables NNGF,GF(t c , At; T 6i Th) and Af^-p G f(^c, At; T 6i Th) factorize: the 
whole time dependence is contained in A|1^ n (t c , At) and is given with a high precision by 
the factor 

, a \ / t c \ sin27rAt/(2T) , . 

f(t c , At) = cos (2^) -^-J^J., (23) 

while the function W(T Gj Th) carries all the information about the energy dependence. De- 
pending on the values of Am 2 from the interval (1) and on T e .Th, the function W(T 6) Th) 
can be positive, negative or zero for E„ = (5.0 — 14.4) MeV; and |W(T e ,Th)| can be greater 
or smaller than 1 (see Table III). Correspondingly, the VO generated and the geometrical 
seasonal asymmetries Aff GF and A|^ s m can have the same or opposite signs and the former 
can be larger or smaller in absolute value than the latter || . The asymmetry which contains 
both the VO induced seasonal effect and the effect due to the geometrical factor, Aq| s , is 
just equal to the sum of the seasonal asymmetry due to the VO only and of the geometrical 
asymmetry, eq. (18). The asymmetry A G p s can be close to zero due to mutual cancelation 
between A^gf an d A^^ ||. This implies, in particular, that one cannot have simulta- 
neously |Aj|gfI ^ Ageom an d IAgFI ^ ^ge^m f° r the same set of values of the parameters. 
However, one of the two asymmetries, A^gf or Aq| s , can be strongly suppressed, while the 
other can have observable values. All these possibilities are realized for the values of Am 2 , 
sin 2 29, E u and T e Th of interest (see further and Figs. 1 and 2). 

Let us note that in ref. |T7]] a Fourier analysis of the predicted time dependence of the 
event rate in, e.g., the Super-Kamiokande detector in the case of the VO solution of the solar 
neutrino problem was performed. The coefficients in the corresponding Fourier expansion 



were found to be expressed in terms of Bessel functions. The authors of ref. [|T7| noticed 
that the constant and the first harmonic term dominate in the expansion. The origin of this 
result becomes clear from eqs. (15) - (18) and (22a) - (22b). 

We expect eqs. (15) - (18) and the above conclusions to be valid for the charged current 
solar neutrino event rates and the corresponding seasonal variation asymmetries to be mea- 
sured with, e.g., the SNO and ICARUS detectors. The expressions for the observables (12) 
and (13) and for the function (19) for the SNO detector can be obtained simply by modifying 
the definition of the "average" in eqs. (5b) and of the averaged quantities in eqs. (10) and 
(19): one has to set i u = and replace in (5b) da(i , e e~)/dT c with d<r(i/ e d — > e~pp)/dT e - 
the differential cross-section of the charged current reaction ^ e + d^e _ +p + p, changing 
as well the lower limit of integration over E u to T + 1.44 MeV. It should also be empha- 
sized that our general results (15) - (19) will not change if one includes in the analysis the 
e~— detection efficiency and energy resolution of a specific detector as long as they do not 
vary with time. 

Let us note that in the case of averaging over a time period not exceeding two months, 
At < T/6, we have 

K(t c , At) = l + 2e cos (2tt^ + 0(~ 1.7 x 10^ 3 ), (24) 
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A— m (t c , At) = 4ecos (2tt|) + 0(~ 3.4 x lO" 3 ), (25) 

and both quantities are practically At— independent. The geometrical asymmetry is maximal 
if t c is chosen to be the time when the Earth is at perihelion: A||^ n (t c = 0, At < T/6) = 4e = 
6.68%. As it follows from eqs. (14) - (18), averaging over a time interval exceeding 3 months, 
At ^ T/4, suppresses the seasonal variation effect. 

3. Enhancing the Seasonal Variation Effect 

We have performed numerical calculations of the observables (12) and (13) for a large 
representative set of values of the parameters Am 2 and sin 2 29 from the VO solution re- 
gion (1) - (2), which are given in Table I. Some of our results are presented in Tables I 
- Ill and in Figs. 1-3. The threshold energy dependence of the event rate ratio (12), 
NNGF(t c , At; T 6j Th), is shown for six selected values of Am 2 and sin 2 29 from the set and 
for t c = (perihelion); T/12; 2T/12; ...; t a = T/2 (aphelion) and At = T/12 in Fig. 1. 
The perihelion-aphelion asymmetries A^ s F GF (t p = 0, At = T/12; T ejTh ) = A^ s F GF (T eiTh ) 
(see eq. (13)) are plotted as functions of T e ,Th in Fig. 2 for all the values of Am 2 and 
sin 2 29 given in the first column of Table I. As Figs. 1 and 2 demonstrate, depending 
on Am 2 and sin 2 29, for T Cj Th — 5 MeV the seasonal effect due to VO can be maximal 
or close to the maximal (Figs, la and 2a, Am 2 = 0.5 x 1CT 10 eV 2 , sin 2 29 = 0.9), or can 
be strongly suppressed (Figs, lb, le, If and 2a, 2c, Am 2 = 0.7; 2.5; 4.4 x 10~ 10 eV 2 ). 
With the increase of T e F h from the value of ~ 5 MeV, the VO induced seasonal effect 
can change drastically. For Am 2 = 0.9 x 10~ 10 eV 2 and sin 2 29 = 1.0 (Figs, lc and 2b), 
for instance, it is negative and close to maximal in absolute value at T e; Th — 5 MeV, goes 
through zero at T e F h — 8.5 MeV and is positive for larger values of T e F h reaching a max- 
imum at T e Th — 12 MeV. If, however, Am 2 = 0.4 x 10~ 10 eV 2 , it decreases monotonically 
with the increase of T e Th (Fig. 2a), while for Am 2 = 1.3 x 10 -10 eV 2 the effect is nega- 
tive and hardly observable at T 0j Th — 5 MeV, but increases in absolute value reaching a 
maximum at T eiTh = 10.9 MeV, A s N c G s F (T CiTh = 10.9 MeV) ^ -14.1% (Figs. Id and 2c). 

In the case of relatively large values of Am 2 = (2.0 — 4.4) x 10~ 10 eV 2 the seasonal effect 
due to VO is relatively small or negligible for 5 MeV < T ejF h & 7.5 MeV (Figs, le, If and 
2c). For Am 2 = 4.4 x 10~ 10 eV 2 and sin 2 29 = 0.9 the VO induced seasonal asymmetry 
increases with T 6j Th reaching a maximum at T 6j Th — 10.1 MeV, A^ GF (T Ci Th = 10.1 MeV) = 
14.8%, then decreases to zero at T Gj Th — 12.2 MeV and becomes negative for larger values 
of T e> Th- The geometrical effect amplifies the asymmetry in the region of the maximum 
and A GF s (T e Th = 10.1 MeV) £ 21.5% (Fig. 2c). At T eiTh = 9.5 (11.0) MeV, we have for 
the indicated values of Am 2 and sin 2 26>: A^ GF = 13.0 (10.0)%. Including the realistic 
experimental conditions in the calculations of the asymmetry will reduce its magnitude, but 
not considerably - by ~ (1-^2)%. The change of T 6j Th from 5 MeV to 10.1 MeV will increase 
the statistical error of the measured value of A G | s (T Gi Th) (see Table I). Nevertheless, values 
of A s N e G s F(GF) (T ejTh ) ~ 15 (21)%, as the one obtained above for Am 2 = 4.4 x 10~ 10 eV 2 at 
T e ,Th — 10.1 MeV, can already be tested with the Super-Kamiokande detector ||. 

In Table II we give the numerical values of A^ GF (T c Th ) at T e Th = 5; 6; 7; 8; 9; 10; 11 MeV 
for the chosen representative set of values of Am 2 and sin 2 29 from the VO solution re- 
gion. For each pair of the latter, the values of A|f GF (T ei Th) corresponding to a maximum 
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of |A^ F (T e , T h)| in the interval T CiTh = (5 - 11) MeV, A s ^ max , and the T CiTh at which 
the maximum is reached, T™f£, are also given. For the set of values of neutrino param- 
eters we have considered, the minimal A^Q S p max is 2.2%; for most of the values, however, 
|A^ G s F max | ^ 4.5%, with the maximal values of the asymmetry reaching ~ (10 15)%. Our 
results show, in particular, that for any given Am 2 from the VO solution region (1) there 
exists at least one value of T e; Th from the interval (5 — 11) MeV, for which the seasonal vari- 
ation effect in the solar neutrino sample of events with recoil electrons having T e > T Cj Th is 
either maximal or very close to maximal (in absolute value). Since the precise value of Am 2 
is still unknown, this implies that one can effectively search for the seasonal effect by forming 
a set of samples of solar neutrino events, which correspond to a sufficiently large number 
of different values of T 6i Th from the above interval, say, for those in Tables I - III, and by 
measuring the seasonal variation in each of these samples. As we have already emphasized, 
the seasonal effect can change dramatically with the sample. 

The correlation between the magnitude of the seasonal asymmetry and the electron 
threshold kinetic energy T e Th discussed above is a manifestation of the fact that the under- 
lying statistical distribution of neutrino events in the case of interest is a bivariate function 
of time and energy, whose marginal distributions are the one-year average electron spectrum 
and the energy integrated time dependence. An example of a possible graphical repre- 
sentation of such a bivariate distribution of the (z/, e~) events is given in Fig. 3a. The 
figure corresponds to Am 2 = 0.7 x 10~ 10 eV 2 and sin 2 29 = 0.85. Each color represents 
a specific value of the one month average event rate ratio Nngf^o At = T/12;T ej Th), eq. 
(12), calculated for t c = (perihelion); T/12; 2T/12; ...; t a = T/2 (aphelion). Other types 
of representations are possible, of course, but the tests we have performed indicate that 
this one is better suited for optimizing the efficiency of the information extraction. These 
maps allow to compare at a glance the predicted time and energy dependence of the ratio 
Nngf^c, At; T Cj Th) for different values of Am 2 and sin 2 29. For example, Fig. 3a clearly 
shows that for the chosen Am 2 and sin 2 29 no observable time modulation is expected at 
T e> xh — 5 MeV, while the maximal modulation takes place for T Cj Th — 9.5 MeV. When 
applied to real data, the presentation technique used to generate Fig. 3a will produce an 
"image" of the bivariate event rate distribution, similar to usual astronomical images ob- 
tained by X-rays or 7-rays telescopes. If the solar neutrino problem is indeed caused by 
vacuum oscillations of solar neutrinos, the event rate distribution "images" will be charac- 
terized by a pattern composed of elongated spots of different intensity, whose distribution 
and shape are functions of the neutrino parameters Am 2 and sin 2 29. It should be possible 
to analyze such images using well tested techniques, as those commonly applied for imaging 
in high energy astronomy in order to ensure not only the presence of relevant patterns, but 
also their consistency with the vacuum oscillations scenario, discarding in this way subtle 
systematical errors induced, for instance, by background time variability. 

As we have shown, for the 8 B neutrino induced event rates of interest and the VO solution 
neutrino parameters from the region (1) - (2) one has: 

N NGF (t c , At; T CjTh ) = q(T C)Th ) + m(T c , Th , At) cos 2tt| + O(10- 3 ) , (26) 

where q is a function of T e Th , m depends both on T c Th and on the sampling interval At, but 
both q and m do not depend on t c . If the data does not exhibit any time dependence, apart 
from that caused by the geometrical effect which is removed from N NGF (t c , At; T e Th ), one 
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should have m(T ej Th, At) = 0. The signature of vacuum oscillations is a significant threshold 
energy dependent deviation of m(T ei Th, At) from zero. For the choice || of the event rate 
normalization, eq. (12), we have q(T eTh ) = 1 for any sampling scheme (see eq. (15)) and, 
as it follows from eqs. (14), (15) and (26), 

m(To , ra , At) . 2£ W(T., Th )?^«, ( 27 ) 

where W(T ej Th) is given in eq. (19). These relations do not necessarily hold true for event 
rate normalizations which differ from that used by us. Equation (26) allows to reduce the 
bivariate distribution N NGF (t c , At; T e!Th ) to the univariate distribution m(T e Th , At) without 
information loss, to the extent that eq. (26) represents a good approximation to the true 
time dependence. In Fig. 3b we compare the time dependence of the one month average 
event rates N NG F(t c , At = T/12;T e Th ) computed utilizing eqs. (5a), (5b) and (12) (points 
with crosses), with those obtained using eq. (26), for T CiTh = 5.0 MeV (solid line), 7.5 MeV 
(dotted line) and 10.0 MeV (dash-dotted line). As Fig. 3b indicates, eq. (26) describes the 
time-dependence of Nngf^c, At = T/12; T Ci Th) with a very high accuracy. This is confirmed 
by a detailed analysis based on the standard Pearson index, performed for Am 2 in the range 
(0.1 -=-4.4) x 10~ 10 eV 2 , sin 2 29 < 1 and T ejTh < 14.0 MeV. The accuracy of eq. (26) increases 
quickly with the increasing of the magnitude of the seasonal asymmetry. The largest errors 
occur for combinations of Am 2 , sin 2 29 and T e xh for which the predicted NGF asymmetry is 
negligible, i.e., is smaller than 0.6%. It is not surprising that eq. (26) is very accurate even 
for Am 2 as large as 4.4 x 10~ 10 eV 2 given the results obtained at the end of the previous 
Section (see eqs. (20) - (22) and the related discussion). 

As we have already indicated, the function W(T e xh) carries all the information about 
the threshold energy dependence of the event rates and the seasonal asymmetries of interest, 
eqs. (15) - (18); its knowledge allows to reconstruct the time dependence of the normalized 
event rates (12). Figure 3c and Table III illustrate the dependence of W on T 6j Th for 
the representative set of values of the neutrino parameters Am 2 and sin 2 29. The function 
W(T e Th ) was calculated numerically by fitting the time-dependence of N NGF (t c , T/12; T e Th ) 
(computed without any approximations) with that of eq. (26), utilizing the least squares 
method. The maximal value of W(T 6i Th) for each pair of Am 2 and sin 2 29, W max , and the 
value of T 6j Th at which it takes place, T™||, are also given in Table III. 

Let us consider an alternative event rate normalization and the corresponding version of 
eq. (26), i.e., Rngf^c, At; T eiTh )/R SSM (T e)Th ) = q'( T c,Th) + m'(T CiTh , At) cos27rt c /T, where 
q'(Te,Th) 7^ 1 represents the "distortion" of the one year average event rate distribution 
and RssM(T ei Th) is the one year average event rate predicted in a standard solar model 
in the absence of neutrino oscillations. Now m'(T ej Th, At) carries information about the 
amplitude of the time variation, while both m'(T ej Th; At) and q'(T e ,Th) determine the en- 
ergy dependence. However, since they are related to the index m(T e) Th, At) in eq. (26) by 
m(T e Th , At) = m'(T e! T h , At)/q / (T e Th , At), it seems more convenient to use the normaliza- 
tion in eq. (12) leading to eq. (26) with only one function of T e)Th (q(T e Th ) = 1), rather 
than having the T 6) Th dependence spread in two functions. 

If the geometrical effect is included in the event rate data, it can be accounted 
for in the observables (12) and (13), as it follows from eqs. (16) and (18), by 
simply adding 1 to the function W(T eTh ) (e.g., in eq. (27)). This is a conse- 
quence of the normalization |J used by us in eq. (12). There exist also simple 
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relations between the function W(T 6j Th) and the different types of seasonal asymme- 
tries considered by other authors. For instance, relation (17) holds for the asym- 
metry A' NGF (At; T CjTh ) = 2(R NGF (0, At; T e , Th ) - R NGF (T/2, At; T eiTh )/(R N GF(0, At; T e , Th ) 
+ Rngf(T/2, At; Te ; Th)- The asymmetry defined in |L8| corresponds to A s = 
2(R GF (t c = 0) - R GF (t c = T/2)/(R GF (t c = T/4) + R GF (t c = 3T/4)) with At = T/4 and 
T c xh = 6.5 MeV and is given by A s = 4e^(l + W(6.5 MeV)). Evidently, eq. (26) can 
be used to relate different seasonal asymmetries and spectral distortion functions discussed 
in the literature, while plots and tables like Fig. 3c and Table III can be utilized to make 
predictions for the latter for a large range of sampling schemes, T 6j Th and values of the 
neutrino oscillation parameters. 

4. Conclusions 

We have studied in detail the threshold energy (T e Th ) dependence of the seasonal vari- 
ation effect in the electron energy integrated solar neutrino signal of the Super-Kamiokande 
detector in the case of the vacuum oscillation solution of the solar neutrino problem. We 
have shown that for the values of the neutrino oscillation parameters Am 2 and sin 2 29 from 
the VO solution region (1) - (2), the time and the T e>Th dependence of the corresponding 
solar neutrino event rate factorize to a high degree of accuracy. This facorization takes a 
particularly simple form, eqs. (15) - (18), when the event rate for a given threshold energy 
T e ,Th is normalized to the one year average event rate for the same threshold energy ||, eq. 
(12). As a consequence of the factorization, the VO induced seasonal variation asymmetry 
is proportional to the geometrical one, the coefficient of proportionality being a function of 
T e xh; but not of time. The asymmetry which contains both the VO induced and the geo- 
metrical seasonal effects, A GF S , is just equal to the sum of the seasonal asymmetry due to the 
VO only, A^ GF , and of the geometrical asymmetry A s ^ m (eq. (18)). The asymmetries A^ G S F 
and Agg^ can mutually cancel and A G | S can be close to zero ||. This implies, in particular, 
that one cannot have simultaneously |An GF | <C and |A G | S | <C A g ^ s m for the same set 

of values of the parameters. However, one of the two asymmetries, A^ GF or A GF S , can be 
strongly suppressed, while the other can have observable values. All these possibilities are 
realized for the values of Am 2 , sin 2 29, E u and T 0j Th of interest. 

The seasonal variation effect exhibits a strong dependence on T eiF h- For the val- 
ues Am 2 = (4.3-4.4) x 1(T 10 eV 2 and sin 2 29 = 0.9, suggested by the current Super- 
Kamiokande data on the e~— spectrum ||, for instance, the VO induced seasonal asym- 
metry is negligible for T Cj Th — (5 — 8) MeV, increases with T ejF h reaching a maximum 
at T CjTh = 10.1 MeV, A s N e G s F (T ejT h = 10.1 MeV) S 14.8%, then decreases to zero at 
T 6i Th — 12.2 MeV and becomes negative for larger values of T ejF h- The geometrical effect 
amplifies the asymmetry in the region of the maximum and A G | s (T ej Th = 10.1 MeV) = 21.5% 
(Fig. 2c). At T C;Th = 9.5 (11.0) MeV, we have for the indicated values of Am 2 and sin 2 29 : 
Af| GF £ 13.0 (10.0)% 5 Results for other VO solution values of Am 2 and sin 2 29 can be 
found in Table II. 



The quoted asymmetry values correspond to ideal experimental conditions; they will be somewhat smaller, 
by ~ (1 — 2)%, if one takes into account the e~ detection efficiency and energy resolution of the Super- 
Kamiokande detector 
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We have shown also that for any given Am 2 and sin 2 29 from the VO solution region (1) - 
(2) there exists at least one value of T e Th from the interval (5 — 11) MeV, for which the VO 
induced seasonal effect in the solar neutrino sample of events with recoil electrons having 
T c > T c Th is either maximal or very close to maximal. Given the fact that the precise values 
of Am 2 and sin 2 26 are still unknown, this suggests that one can effectively search for the 
seasonal effect by forming a set of samples of solar neutrino events, which correspond to a 
sufficiently large number of different values of T 6j Th from the above interval, say, for those 
in Tables I - III, T cTh = 5; 6; 7;...; 11 MeV, and by measuring the seasonal variation in 
each of these samples. The seasonal effect can change dramatically with the sample. The 
results summarized above are illustrated in Figs. 1-3 and Tables II - III. Although they 
have been derived for the Super-Kamiokande detector, we expect similar results to be valid 
for the SNO and ICARUS detectors. 
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TABLES 



TABLE I. The ratio R(T Ci Th) /R(5 MeV) of the one year average event rates due to 8 B neutrinos 
in v — e~ detectors (Super-Kamiokande, etc.) for several values of T ei Th and for a set of values of 
Am 2 and sin 2 26 from the VO solution region. 
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TABLE II. The perihelion - aphelion asymmetry Afgf F (t c = 0, At = T/12; T CjTh ) 
= Af|Q F (T ej Th) in % for different values of T Ci Th and a set of values of Am 2 and sin 2 29 from 
the VO solution region. Given also are the values of the asymmetry corresponding to a maximum 
of |A|Jq s f | in the interval T Cj Th = (5 — 11) MeV, A^ s F max , and the T 6) Th at which it is reached, 
TJ£f£. The asymmetry including the geometrical effect AgfP (0, T/12; T C)Th ) = Af|g F (T CjTh ) + 4e. 
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4.43 


-0.68 


0.86 


2.33 


3.43 


4.08 


4.38 


4.42 




ooc 

bo bo be 
ooc 


0.65 
0.70 
0.85 


11.00 
11.00 
11.00 


2.82 
3.28 
5.19 


-1.38 
-1.60 
-2.56 


-0.47 
-0.55 
-0.91 


0.56 
0.66 
1.11 


1.46 
1.72 
2.92 


2.14 
2.51 
4.18 


2.58 
3.01 
4.90 


2.82 
3.28 
5.19 




0.85 
0.85 


0.75 
1.00 


11.00 
11.00 


3.54 
8.99 


'-! in 
-6.73 


-2.06 
-5.18 


-0.64 
-1.78 


0.84 
2.46 


2.08 
6.00 


z.yy 
8.15 


'-! 

O.04 

8.99 




0.90 
0.90 


0.75 
1.00 


5.00 
5.29 


-3.89 
-8.22 


-3.89 
-8 12 


-3.33 
-7.85 


-2.11 

-5.64 


-0.58 
-1.70 


0.90 
2.74 


2.13 
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2.99 

8 35 




1.00 
1.00 


0.75 
1.00 


5.86 
6.72 


r — i \ r 7 

-5.07 
-10.56 


-4.75 
-8.68 


-5.06 
-10.17 


-4.61 
-10.49 


-3.49 
-9.13 


-1.99 

-5.85 


-0.39 
-1.22 


1.03 
3.29 




1.10 


1.00 


8.25 


-12.24 


-7.61 


-9.77 


-11.45 


-12.23 - 


11.66 


-9.33 


-5.36 




1.20 


1.00 


9.46 


TTJ — T7T 

-13.40 


-5.88 


-8.18 


-10.39 


-12.15 - 


13.21 


-13.24 


-11.83 




1.30 


1.00 


10.85 


— i — i — i — p 
-14.11 


-3.98 


-6.15 


-8.53 


-10.69 - 


12.45 


-13.68 


-14.09 
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0.90 


9.79 
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14.31 


0.10 


-0.12 


0.87 


7.10 


13.18 
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4.40 


0.90 
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14.75 


0.92 


1.04 


2.34 


1.62 


8.59 


14.69 
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TABLE III. Values of the 


function 


W(T e , Th 


) for various T c 


Th and A 


m 2 and 


sin 2 29. 




(10 


Am 2 
-io eV2) 


sin 2 29 


rpmax 
L e,Th 

(MeV) 


w max 


5.0 


6.0 


T 

7.0 


3 ,Th (MeV) 
8.0 9.0 


10.0 


11.0 




0.40 


1.00 


5.00 


0.61 


0.61 


0.56 


0.51 


0.46 


0.41 


0.37 


0.33 




0.50 
0.50 


0.60 
0.90 


6.17 
5.51 


0.33 
0.68 


0.32 
0.67 


0.33 
0.67 


0.32 
0.64 


0.31 
0.60 


0.29 
0.55 


0.27 
0.50 


0.25 
0.45 




0.55 


1.00 


6.17 


0.97 


0.91 


0.97 


0.95 


0.88 


0.81 


0.73 


0.66 




0.60 
0.60 


0.70 
0.90 


8.02 
7.30 


0.47 
0.80 


0.34 
0.61 


0.42 
0.74 


0.46 
0.79 


0.47 
0.79 


0.46 
0.75 


0.44 
0.70 


0.42 
0.64 




0.65 


0.85 


8.49 


0.73 


0.38 


0.56 


0.68 


0.73 


0.73 


0.Y0 


0.66 




0.70 


0.85 


9.50 


0.76 


U. 10 


0.40 


0.60 


0.71 


0.75 


U. / 


11 7'-! 
U. 1 o 




0.75 


0.80 


10.63 


0.67 


-0.10 


0.13 


0.36 


0.52 


0.62 


0.67 


0.67 




ooc 

bo bo be 
ooc 


0.65 
0.70 
0.85 


11.00 
11.00 
11.00 


0.43 
0.50 
0.79 


-0.21 
-0.24 
-0.39 


-0.07 
-0.08 
-0.13 


0.09 
0.10 
0.17 


0.22 
0.26 
0.45 


0.33 
0.38 
0.64 


0.39 
0.46 
0.75 


0.43 
0.50 
0.79 




0.85 
0.85 


0.75 
1.00 


11.00 
11.00 


0.54 
1.37 


-0.45 
-1.02 


-0.31 
-0.78 


-0.09 
-0.26 


0.13 
0.38 


0.32 
0.92 


0.46 
1.25 


0.54 
1.37 




0.90 
0.90 


0.75 
1.00 


5.02 
5.41 


-0.59 
-1.24 


-0.59 
-1.23 


-0.50 
-1.19 


-0.32 
-0.85 


-0.08 
-0.24 


0.14 

0.43 


0.33 
0.97 


0.46 
1.28 




1.00 
1.00 


0.75 
1.00 


5.94 
6.78 


-0.77 
-1.60 


-0.72 
-1.32 


-0.77 
-1.54 


-0.70 
-1.59 


-0.53 
-1.38 


-0.30 
-0.87 


-0.05 
-0.17 


0.16 
0.52 




1.10 


1.00 


8.17 


-1.85 


-1.15 


-1.48 


-1.73 


-1.85 


-1.76 


-1.41 


-0.79 




1.20 


1.00 


9.50 


-2.03 


-0.90 


-1.24 


-1.58 


-1.84 


-2.00 


-2.00 


-1.78 




1.30 


1.00 


10.83 


-2.14 


-0.61 


-0.94 


-1.30 


-1.62 


-1.89 


-2.07 


-2.13 




2.50 


0.90 


9.68 


2.19 


0.02 


-0.02 


0.14 


1.10 


2.04 


2.14 


1.73 




4.40 


0.90 


10.05 


2.30 


0.14 


0.16 


0.37 


0.25 


1.35 


2.29 


1.50 
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FIGURE CAPTIONS 



Figure 1. The t c and T C Th dependence of the normalized one month average event rate 
Nngf^o At = T/12;T ej Th)> eq. (12), in v — e~ detectors (Super-Kamiokande, etc.) for six 
representative values of Am 2 and sin 2 29 from the VO solution region, given in Table 1 
under the numbers: n = 3 (a), 8 (b), 16 (c), 21 (d), 22 (e), 23 (f). The solid, dotted, dash- 
dotted, dashed lines and the lines with squares, diamonds, stars in each of the sub-figures 
(a) - (f) correspond to seven equally spaced values of t c = (perihelion); T/12; T/6; T/4; 
T/3; 5T/12; T/2, respectively. The seasonal asymmetry is given by Aj|Q F (t c , T/12; T 6i Th) 
= N NGF (t c ,T/12;T CiTh ) - N NGF (T/2-t c ,T/12;T CiT h), < t c < T/4, and is maximal for 
t c = 0. 

Figure 2. The dependence of the perihelion - aphelion asymmetry Aff GF (t c = 0, At = T/12; 
T c ,Th) (in %), eq. (13), on T 6j Th for different values of Am 2 and sin 2 29 from the VO solution 
region. The curves in the panels (a), (b) and (c) correspond to the values of Am 2 and sin 2 29 
in Table I numbered, respectively: 1 - 8, 9 - 16 and 17 - 23. Each curve of a given type (solid, 
dashed, dash-dotted, etc.) is obtained for the values of Am 2 and sin 2 29 which are marked 
with the symbol of the curve (solid line, dashed line, dash-dotted line etc.) in the second 
column of Table 1. The asymmetry which includes the geometrical effect Aq F s (0, T/12; T e F h) 
= A S N C G S F (0, T/12; T CjTh ) + 4e, 4e = 6.68%, is also shown. 

Figure 3. a). Am "image" of the one month average event rate ratio (dis- 
tribution) N NG F(t c , T/12;T 6i Th) (eq. (12)). Each vertical stripe corresponds to 
t c = 0; T/12; T/6; T/4;...; T, as is indicated on the horizontal axis. The vertical scale 
on the right-hand side allows to convert a color into a value of the ratio, b) The event rate 
ratio shown in a), plotted as a function of cos27rt c /T for T CiTh = 5.0 MeV (solid line), 7.5 
MeV (dotted line), 10.0 MeV (dashed-dotted line). The lines are the best fits to the com- 
puted ratio indicated by "+". c) The dependence of W(T ejTh ) on T CjTh . All results shown 
in the figure are derived for Am 2 = 0.70 x 10~ 10 eV 2 and sin 2 29 = 0.85. 
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